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Phonon decay in silicon nanocrystals 

A. A. ProkoficvQ A.N. Poddubny, and I.N. Yassievich 

lojfe Physical- Technical Institute of RAS, 194021 St. Petersburg, Russia 

(Dated; April 8, 2013) 

The decay of the optical phonon into the two phonons of smaller energy is calculated for Si 
nanocrystals. The rate of the process is in the range of 1 to 10 ps. Such anharmonic phonon decay 
may control the energy relaxation rate of excited carriers in Si/Si02 nanocrystals. Relevance of the 
phonon decay to the experimentally observed hot carrier photoluminescence is discussed. 



during 1 -^ 100 ps after the pump pulsei^ii To resolve the 
paradox we propose the following seenario, see Fig. [1] (i) 
Eleetrons and holes emit optieal phonons very fast, (ii) 
These optical phonons do not escape the nanocrystal due 
to their low group velocity and differences between the 
phonon spectra in Si and Si02 .— Instead, they are quickly 
reabsorbed by the confined carriers, so that the total en- 
ergy of the eleetron-phonon system is conserved, (iii) 
Phonon recycling process stops after the only after the 
anharmonic decay of high-energy phonons into the lower 
energy ones. This low energy acoustic-like phonons can 
relatively easy leave the nanocrystal, so the energy of the 
system "carriers-|-phonons" is dissipated. Thus, the en- 
ergy relaxation rate is governed not by the phonon emis- 
sion rate, but by the optical phonon decay rate. Similar 
mechanisms arc known in literature and important, for 
instance, for GaN devices i^ 

The paper is organized as follows. Calculation ap- 
proach for the phonon modes is outlined in Sec. |lll 
Phonon emission and phonon decay processes are con- 
sidered in Sec. IIIII and Sec. IIVI Main paper results are 
summarized in Sec. [V] 



II. PHONON MODES OF THE SI 
NANOCRYSTAL 

Dispersion of phonons in bulk silicon can been de- 
scribed within various Keating-type modelS)^"— where 
the atoms arc effectively treated as balls connected by 
springs. In our approach we adopt the simplest two- 
parameter Keating model from Ref. Il2l which allows 
straightforward generalization for anharmonicity. Ac- 
cording to this model, the strain energy is given by 




FIG. 1. The schematical representation of the suggested re- 
laxation mechanism. Confined carriers in Si nanocrystal emit 
optical phonons (as well as absorb them) very fast. Opti- 
cal phonons decay into acoustic ones, which can leave the 
nanocrystal taking the energy away. 



I. INTRODUCTION 

Silicon nanocrystals form a rapidly developing area of 
research with a number of promising applications for op- 
toelectronics and photovoltaics. Rapt attention is now 
attracted to such processes as space-separated quantum 
cuttingir— and hot photoluminescence;^ which intrinsi- 
cally involve hot confined carriers, optically pumped in- 
side the nanocrystal. This raises the demand to clarify 
the mechanisms, controlling the hot carrier lifetime. 

The general understanding is that the energy relax- 
ation of hot carriers is determined by the emission of 
optical phononsi^i^ For lowest excited states this process 
can be slowed down in nanocrystal with respect to that 
in the bulk due to the energy spectrum discretization.— 
However, this argument does not hold for hot excited 
states, which require completely different consideration. 
We will show below, using the atomistic calculations of 
carrier and phonon spectra, that the rate of energy re- 
laxation due to optical phonon emission may reach the 
values up to 0.1 eV/fs. Such ultrafast relaxation is in 
clear contradiction to the various experimental manifes- 
tations of the hot carriers presence in the nanocrystal 
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where oq is the lattice constant, and A denotes the 
change to the unstrained lattice: A(r^) = r^ — Tq, 
l^ivij • r,fc) = (r^j -Yik) - {r°j • r°j.). The two constants 
a and (3 describe the bond-stretching and bond-bending 
forces. In the harmonic approximation the Hamiltonian 
of the vibrational system assumes the form 
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FIG. 2. Calculated density of phonon states in Si nanocrystal 
with diameter 2.5 nm. Open bars correspond to "optical" 
phonons (with energies above the phonon energy gap). The 
inset shows phonon dispersion of bulk silicon calculated in the 
same model. 



where Mgi is the mass of silicon atom, p is the momentmn 
operator, and 



V,., 
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(3) 



The phonon modes in the nanocrystal have been calcu- 
lated using the clamped boundary conditions: the atoms 
displacements of all the atoms outside the nanocrystal 
were set to zero. Relaxation of the nanocrystal lattice 
has been neglected for simplicity. The values of con- 
stants a and /? from Ref. [l2 have been used. The cal- 
culated density of phonon states in Si nanocrystal with 
diameter 2.5 nm is shown in Fig. [2l The figure demon- 
strates that the phonon modes can be divided into two 
distinct groups. First group corresponds to the acoustic- 
like modes with the energy hw < 40 meV (filled bars). 
Size quantization of the phonons explains the absence 
of the modes with the energies less than 12 meV. Sec- 
ond group of the modes with fiw > 40 meV corresponds 
to optical-like phonons and is manifested as two distinct 
peaks in the density of phonon states (outlined bars). 
These results are in qualitative agreement with Ref. Il3| 
where a detailed theoretical study of phonon modes in Si 
nanocrystals is presented. 



III. PHONON EMISSION PROCESSES 

Phonon emission processes, usually responsible for en- 
ergy relaxation, can be considered as follows. First, the 
energy spectra of confined carriers and phonons in a 
nanocrystal are calculated separately. Second, the in- 
teraction of carriers with phonons is analyzed using the 
fact that the distortion of the atomic positions due to 
the vibrations changes the matrix elements of the tight- 
binding Hamiltonian for the confined carriers. 



We have used nearest-neighbor sp^d^s* tight-binding 
method to calculate the probability of single-phonon 
transitions of excited electrons and holes in Si 
nanocrystalsJ^ The dangling bonds at the nanocrystal 
surface were "passivated" with hydrogen atoms, the pa- 
rameters were taken from Ref. |l^ The relaxation of the 
nanocrystal lattice from the bulk has been neglected. 
Calculated densities of states of confined carriers and 
more details of the model can be found in Ref. \1Q. 

The emission rate of the single optical phonon has been 
calculated according to the Fermi Golden rule as 



7 = T E E</l^i-PhN)''5(i?, -Ej- tu.^) , 



(4) 



where Ei and Ef are the energies of initial and final state 
and the summation runs over the phonon states with the 
energies fuou. In numerical calculation the 5- function, 
accounting for energy conservation, has replaced by a 
Lorentzian with HWHM equal to 1 meV. The Hamilto- 
nian of carrier interaction iJc-ph with the given phonon 
mode v has been found by taking the derivative of the 
tight-binding Hamiltonian over displacement: 
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where 5r\" = uY' — u".^' are the atomic displacements, 
corresponding to the phonon mode v. The dependence 
of the matrix elements of the tight-binding Hamiltonian 
on the interatomic distance has been taken from Ref.—. 
Energy relaxation rate for a given state i is obtained by 
multiplying each term in Ec^. ^ by hw,^. 

Fig. [3] shows calculated energy relaxation rates of elec- 
trons and holes in 2.5 nm silicon nanocrystal due to 
single-phonon transitions from all possible excited states. 
The figure demonstrates that the energy relaxation rate 
can be quite fast and may reach the values ^ 10^^ eV/s, 
i.e. the energy of one eV could be lost within the tens 
of femtoseconds. Majority of the transitions correspond 
to a somewhat slower relaxation rate within the range 
~ lO^^-f-lO^"^ cV/s. These new relaxation rates are signif- 
icantly higher than those reported in our previous workj^ 
the discrepancy is due to the technical mistake in Ref. [la. 

Noteworthy, the fastest relaxation processes are due 
to the emission of optical phonons, see open symbols in 
Fig. [31 Thus, hot carriers in nanocrystal may generate 
a large population of optical phonons. If these phonons 
do not leave the nanocrystal fast enough, they can be 
reabsorbed by the carriers. This means that the phonon 
emission process itself does not lead to the dissipation of 
the total energy in the system. The energy relaxation 
is controlled by another processes, such as anharmonic 
phonon decay of optical phonons into acoustic ones, an- 
alyzed below j^iiiii^ 
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FIG. 3. Relaxation rate of electrons and holes (negative en- 
ergies) in Sises nanocrystal {d « 2.5 nm). Each point corre- 
sponds to single-phonon transition rate from the state with 
given energy, related to the energy of phonon. Open sym- 
bols indicate the transitions involving phonons with energy 
higher than 40 meV. The histogram in panel (b) referres to 
the density of electronic states in the nanocrystal. 



IV. PHONON DECAY PROCESSES 

Phonon anharmonicity may be accounted for by 
adding the terms in the strain energy Eq. ([T]), cubic in 
the atomic displacement. This is done by including the 
dependence of the force constants a and /? in Eq. ([J) on 
the bond lengths rjk<^ 



ao 



/5 = /5o 






'_jk_ 



1 12 




111 



(6) 

(7) 



Here, the values r^^, correspond to the undistortcd lattice, 
and ao,/3o are the values of constants a, /3 used in har- 
monic approximation. The anharmonic part of the strain 
energy, cubic in displacement, can be then reduced to 
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Using the second-quantization formalism, we rewrite the 
harmonic part of the phonon Hamiltonian as 



H = J2t^u(b% + l 



(9) 



The terms of the anharmonic Hamiltonian responsible for 
the decay of the phonon mode v into the modes fi and k 
reads 



^anh — 7 ^ -^vnKMiib^bi^ , 
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FIG. 4. Phonon decay rate in Siaea nanocrystal (d « 2.5 nm). 
The density of phonon states in the nanocrystal is shown in 
panel (b) for reference. 



where 
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The decay rate of the phonon mode a reads 



(11) 



(12) 



In numerical calculation the ^-function, accounting for 
energy conservation, has been replaced by a Lorentzian 
with HWHM equal to 0.15 meV. 

The calculated rates of phonon decay for all the modes 
of Sises nanocrystal with the diameter close to 2.5 nm 
are shown in Fig. 21 The figure demostrates that the 
decay time monotonously decreases with the phonon en- 
ergy and reach the value ^ 1 ps. This growth of the 
phonon decay rate with phonon energy reflects the in- 
creasing density of final states. Our analysis indicates, 
that the phonon emission rate weakly depends on the 
nanocrystal size. It may be estimated using the following 
simple argument. The element of cubic phonon-phonon 
interaction is approximately given by the product of the 
phonon energy and the ratio of the vibration amplitude 
to the lattice constant. 
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The value of u^, can be found from the normalization con- 
dition AfgiM^ = huj^. Finally, the density of final phonon 
states in Eq. (|12p can be estimated as 1/wi/. This yields 
(t^)~i w 10^^ s~^, in agreement with Fig. |4l 



V. CONCLUSION 

Theory of phonon emission and phonon decay pro- 
cesses in silicon nanocrystals has been developed. The 



calculation has been carried out using the empirical 
sp^(fs* model for the states of confined carriers and elec- 
trons combined with the anharmonic Keating model for 
the vibrations of the nanocrystal atoms. 

The calculations yield very short times of optical 
phonon emission, reaching tens of femtoseconds for high 
energy optical-like phonons. This is an indication for 
a strong clectron-phonon interaction. Hence, emitted 
phonons can be quickly reabsorbed, and the total energy 
of the system "electrons-|-phonons" is conserved. The en- 
ergy dissipation rate is then controlled by the anharmonic 
decay of the high-energy phonons into several low-energy 
acoustic-like ones, which then leave the nanocrystal. 



Calculated phonon lifetimes due to three-phonon decay 
processes are on the order of tens of picoseconds and 
weakly depend on the nanocrystal size. Such decay can 
govern the total energy relaxation rate of hot carriers in 
silicon nanocrystal. The value of the rate corresponds to 
the lifetime of hot carriers determined from the decay of 
the hot photolumincsccncc in Si/Si02 nanocrystalsi^ 
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